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SUMMARY 

This paper investigates the feasibility of a proposed technique for detecting delaminations using piezo- 
electric layers or patches embedded or bonded to a composite structure. Variations in the voltage generated in 
the piezoelectric layers indicates the presence and location of delamination, when the structure is excited either 
externally or via piezoelectric actuators. The theoretical foundations of a method for predicting the dynamic 
response of delaminated composite beams with piezoelectric layers are described. The governing equations are 
presented for the case of external vibroacoustic excitation, as well as, for the case of locally induced vibrations 
by some of the embedded piezoelectric elements. An exact solution is developed within the limits of linear 
laminate theory. Applications illustrate the feasibility of delamination detection in cantilever beams. The results 
illustrate that the proposed technique may provide accurate detection of the presence, size and location of a 
delamination. 


INTRODUCTION 

Since the advent of composite materials, the development of delamination cracks have represented one of 
the principal concerns to their engineering applications. Delamination cracks may result in catastrophic failure, 
but also may easily remain undetected. Consequently, numerous theoretical and experimental studies have been 
conducted to elucidate the effects of delaminations on the static and dynamic response of composite structures. 
Additional research has addressed the nondestructive detection of delaminations. 

The problem of delaminations has received considerable attention as a fracture mechanics problem and the 
recent papers of Chatteijee (ref. 1), Pindera (ref. 2), Nilsson and Storakers (ref. 3) and Thangjitham and Choi 
(ref. 4) are mentioned. The vibration and buckling of delaminated composite structures has been also investi- 
gated. Simitses (ref. 5) provided a comprehensive bibliography on the subject of buckling of composite plates 
and beams and a typical formulation of a delaminated beam problem based on a subdivision of the beam into 
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four regions. Anastasiadis and Simitses (ref. 6) reported also analytical treatments for the case of weak bonding 
between the delaminated plies. Lee et al. (ref. 7) applied a layer- wise theory to the buckling of a composite 
beam and considered multiple delaminations through the thickness of the beam using a finite element 
formulation. 

A large number of nondestructive, vibration based, techniques for detection of delaminations have been also 
developed. A review of vibration methods used for nondestructive evaluations (NDE) of laminated structures 
was published by Cawley and Adams (ref. 8). Cawley (ref. 9) also reviewed a number of methods suitable for 
low-frequency NDE. Tracy and Pardoen (ref. 10) presented theoretical and experimental vibration studies on the 
effects of delamination on the natural frequencies of specialty composite beams. They subdivided the beam in 
four regions, similar to the approach employed in the buckling studies, used traditional Euler beam theory and 
neglected extension-flexure coupling effects. Effects of the delamination on the higher natural frequencies were 
observed, depending on the location and size of delamination. Nagesh Babu and Hanagud (ref. 11) and Paolozzi 
and Peroni (ref. 12) used finite element models to study the same problem. Their results confirmed the 
conclusions of Tracy and Pardoen (1989) regarding a reduction in natural frequencies as a result of delamina- 
tion. Additional finite element solutions have been published by Tenek et al. (ref. 13) (1993). They concluded 
that the effect of delaminations increases for higher natural frequencies. Shen and Grady (ref. 14) developed 
equations of motion for a delaminated beam using the Hu-Washizu variational principle and Timoshenko beam 
theory, and also reported experimental results and correlations. Saravanos (ref. 15) developed an exact analytical 
procedure for the analysis of delamination effects on the natural frequencies, mode shapes and modal damping 
of composite beams. The last approach, that represents a version of a layer-wise theory, allows tracing a number 
of delaminations through the thickness of the laminate. 

“High-frequency” nondestructive methods have been also proposed for delamination detection and are 
briefly reviewed here for the shake of completeness. Velocity measurements of surface waves have been used to 
nondestructively detect various defects, including delaminations (ref. 16). Buynak et al. (ref. 17) used an 
ultrasonic technique to produce images of delaminations in graphite-epoxy laminates. Acousto-graphic imaging 
techniques (ref. 18), and a photothermal method was discussed (ref. 19). Acoustic emission measurements were 
proposed (refs. 20 and 21) to monitor propagation of delaminations. Thermo-acoustic emission was considered 
(ref. 22) and shearography techniques have been also proposed (ref. 23). Moire interferometry (ref. 24) was also 
used to detect strain concentrations associated with the presence of delaminations in graphite-epoxy laminates. 

A relatively new area of non-destructive delamination monitoring is associated with adaptive or smart 
structures. These structures may incorporate piezoelectric materials, shape memory alloys, electro-rheological 
fluids, electrostrictive, and magnetostrictive materials. Nagesh Babu (ref. 1 1) and Hanagud et al. (ref. 25) 
recommended the use of piezoelectric sensors to detect delaminations in composite beams. An analytical model 
for automatic detection and control of delaminations growth using piezoelectric sensors and actuators was 
discussed (ref. 26), and a four-region model of a delaminated beam based on Timoshenko beam theory was 
used, while solutions were obtained by the finite element method. Teboub and Hajela (ref. 27) also implemented 
a four region model of a delaminated beam with piezoelectric sensors for delamination detection under static 
loads using Kirchhoff type assumptions in each region. An experimental and numerical investigation on 
applications of piezoelectric sensors to detect delaminations in composite laminates was carried out (ref. 28). 

Kim et al. (ref. 29) proposed the use of strain gages bonded to the surface or embedded within the laminate. 

The present paper investigates the feasibility of a proposed detection technique based on monitoring the 
spatial voltage distributions across distributed piezoelectric sensors. The proposed technique may provide simple, 
real-time delamination detection capabilities and requires minimal instrumentation. Hence, it may enhance the 
health monitoring capacity of composite structures and improve their reliability. Mechanics for composite 
piezoelectric beams with multiple delaminations and piezoelectric patches bonded to the top and bottom surfaces 
are developed. An analytical solution is presented, within the limits of the assumptions of linear theory of elastic 


2 



composite beams. The solution is further concentrated on the application of piezoelectric sensors to detect 
delaminations. Evaluations of the method on cantilever beams illustrate the merits and effectiveness of the 
proposed technique. 


Theoretical Background 

This section presents the theoretical formulation for modelling the dynamic response of delaminated 
composite beams with embedded piezoelectric sensors and actuators. A number of issues can be addressed in 
the framework of this formulation including, the detection of delaminations and the control of the dynamic or 
static response of delaminated beams using some of the piezoelectric elements as sensors and the remaining 
piezoelectric elements as actuators. 

Consider a multilayered beam formed from composite plies and piezoelectric layers, as shown in figure 
1(a), with delamination cracks between the composite plies. The delaminations may be at any axial location and 
may overlap each other. The piezoelectric patches may be bonded to the surfaces or embedded between the 
composite plies of the beam (fig. 1(a)). Although the theory can handle embedded layers, additional consider- 
ation may be required in this case due to local stresses in the vicinity of piezoelectric patches (refs. 30 to 32). 

For analytical purposes, the beam may be subdivided into four typical segments. A, B, C, and D respective- 
ly, as shown in figure 1(b). Type A corresponds to an “ordinary” beam without de laminations and patches. Seg- 
ments B correspond to the portions of the beam with piezoelectric patches but without delaminations. The 
segments of the beam with delaminations but without patches are denoted as type C. Finally, the segments D 
include both delaminations and piezoelectric layers. Note that because of the delamination, the segments C and 
D include two additional subregions separated by a delamination. According to the approach used in this paper, 
equations of motion are formulated for each subregion within the segments C and D and for each segment 
(region) A and B independendy. The coupling between the admissible solutions is introduced by the continuity 
conditions for displacements, forces and moments at the boundaries of each region. 


The three-dimensional constitutive relations for a piezoceramic material with anisotropic in-plane elastic 
properties and a macroscopic polar axis aligned with the z-axis are (ref. 33): 
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where e and o indicate engineering strains and stresses respectively. The matrices s E mn , d mn and , T are the 

nm ran m 

elastic compliances, piezoelectric constants, and dielectric permittivity respectively. Superscripts E and T 
indicate that the corresponding properties should be measured in constant electric field and constant stress 
conditions respectively. E i and Dj are the components of the electric field and electric displacements. In the case 
of a thin beam it may be assumed that c y = o z = 0^= 0, which simplifies the constitutive equations as follows, 
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where c* n = 1/S U is the equivalent axial modulus of the material. 


It is assumed that the problem is geometrically linear, moreover, the transverse shear deformations and 
rotary inertia are considered negligible. It is further assumed that u k = u k (x) and w k = w k (x) are the axial and 
through-the-thickness displacements at the mid-plane of the k-th region (or subregion), while the transverse 
displacement is assumed negligible (v^x)^). Then the strain-displacement relationship for thin sections is: 
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where z k is the distance from the mid-plane axis of each region. 


Integration of the previous equation through the thickness of the beam and combination with equations (2) 
yields the stress and moment resultants N x k and M x k respectively, in the k-th region of the beam. 
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where A, B, and D are the extensional, coupling and bending stiffness of the beam, 
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and F, G are the piezoelectric force and moment resultants respectively, 
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n lk ^ 

m- 1 , 

‘-m 

indicates the number of plies across the k-th region (or subregion) of the beam, In segment A (see fig. 1(b)), 
n lk is the number of composite plies through the depth h of the beam; in segment(s) B n^ must be modified to 
incorporate the piezoelectric patches; in segments C and D, the limits of integration must be chosen to include 
the corresponding delaminated sublaminates. 

The equations of motion in the k-th region of the beam, assuming that the axial inertia effects are 
negligible, are: 
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where q(x,t) is an external load, and 
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where p is the linear density of the beam. 

Substitution of equations (4) into equations (7) yields the equation of motion for the kth region: 

A k u X' B k w X*~ F k( E z^ '■ 0 


( 9 ) 
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( 10 ) 


Equations (9) and (10) can be used to analyze the motion of the beam when it is excited by external loads 
(sensory beam), as well as, when it is excited by some of the piezoelectric layers (sensory/active beam). In the 
sensory case, an external load q(x,t) is present, while the imposed electric field E z in the piezoelectric layers is 
zero. In a sensory/active beam, electric fields E z may be imposed in some of the piezoelectric layers to induce 
additional vibrations, while external mechanical loads q(x,t) may be present. 

The solution must satisfy the boundary conditions at the supports. For example, the free-end conditions are 
N x = M x = V x = 0, where V x is the shear force, and the clamped conditions are u= w= w x = 0. Additional 
kinematic continuity conditions between the k-th and (k+l)-th regions at x = Xj. are imposed: 

W k \ X r) ' w * 'Kl’ (11) 
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Continuity in the balance of axial forces, bending moments and shear forces is also required 
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where the indices k and j identify the number of subregions to the left (x' r ) and to the right (x + r ), respectively, 
of a cross-section at x = x r The coordinates z k reflect the distance of the elastic axes of the corresponding 
subregions to the mid-plane of the beam. The shear forces in equation (12.3) are given by: 
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where (r=k or j). In the absence of external axial forces, the last term in the right side of equation (13) may be 
neglected and equation (12.3) becomes: 

E M = E M i*W (I 4 ) 

k j 

The governing equations (9-10) together with the boundary conditions and continuity conditions (12) 
formulated in this section represent the necessary analytical relations for the solution of problems of active and 
passive control of composite beams with delaminations. In the next section applications of the theory for the 
detection of delaminations are presented. 


Detection of Delaminations 

This section further extends the method for the detection of delaminations under dynamic excitation. In the 
proposed scheme, distributed piezoelectric sensors are used to detect the presence, size and location of 
delaminations, when the beam is excited either passively or actively by some of the piezoelectric layers. In both 
cases, the voltages generated across the piezoelectric sensors in an intact (pristine) beam, are different than those 
in the delaminated counterpart. 

Consider the more general case where the beam is excited by either an external harmonic load g^ x y tot or 

an applied harmonic electric field E z ( x)e mt at some of the piezoelectric layers (actuators). It is customary to 

neglect the effect of a self-generated voltage on the motion (ref. 34), however, techniques for incorporating 
these effects have been reported by Heyliger and Saravanos (ref. 32) and may be included in future work. The 
displacements are assumed to have the form 

(«*, w k ) - \u k ,W k \e sx e iw (15) 

and the equations of motion (9-10) are reduced to: 
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where A^is a constant and q k is the amplitude of the external load acting on the region. 
The solution of equations (16) is: 

nr k s' k • S' k S' k S lc X , s' k 

W - C j smyi C 2 cos s k x+ C 3 e + C 4 


(16.2) 


J- k ^k 

f/ *= c/+ _ix* 


pit , 


where C t k through C 5 k are coefficients of integration and 


s k z 


Pfc® 


P, - 


e Q k >R k 

(17.1) 

3 V**- c 4 V 5 *V p k 

(17.2) 


(18.1) 

x )/p*<o 2 

(18.2) 


(18.3) 


D k - B;/A k 

ijf t (£ 2 )*, R t - G t (_ w 

e» ■ - 

Note that integration of the system of equations (16) results in 5 constants of integration for each beam region 
including the constant N k . The number of boundary and continuity conditions is proved to be always equal to 
the number of coefficients of integration. Hence, a linear system of algebraic equations results by substituting 
equation (17) into the boundary and continuity conditions, of the form 

W M ■ W (19) 

where [a ;j ] is a square matrix, [G] is the vector of unknown coefficients in equations (17), and [fj] is a vector 
of loading terms. The elements or the matrix [a] are independent to the method of excitation (passive or active). 
Once the constants of integration are determined, the displacements, strains and stresses can be evaluated. 


The electric potential difference (voltage) 5<t> generated across the sensory piezoelectric layers or patches 
can be calculated as 

8 <K0 ■ - J E z {i)dz (20) 

h P 

where hp is a thickness of the corresponding piezoelectric layer. The positive direction of the electric field 
component coincides with the one of the z-axis. Combining equations (2), the electric field can be obtained 
from 
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In a sensory piezoelectric layer the electric circuit remains practically open, hence, the total charge D z over 
the length Ip of the electrodes remains constant and can be assumed equal to zero, 


f D. dx = 0 (22) 

Ji p 

Combining equations (21) and (22), and considering that the terminals impose a constant electric potential 
condition along its length 1 , the voltage output across the electrodes of a sensor can be obtained as: 
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Passive Excitation .-In the case of external mechanical excitation, the terms in equation (18.2) vanish, hence, 
the right hand side term {fj} in equation (19) contains only terms related to the mechanical load. 

Piezoelectric Excitation - In addition to external mechanical excitation, voltage may be applied across some 
of the piezoelectric layers/patches to actively excite the beam. This inclusion of piezoelectric excitation requires 
only modification of the loading vector {fj} in equation (19). If the motion of the beam is excited by a uniform 
sinusoidal voltage applied to a continuous actuator-layer, the solution of the equations of motion is given by 
equations (17) (where F k may or may not be zero) and U k is complemented by the following additional term: 

p k - £; c,\4dz < M > 

k z* 

APPLICATIONS AND DISCUSSION 
Materials and Assumptions 

Evaluations of the previously described model were performed for sensory composite beams. The composite 
plies were T300/934 graphite/epoxy of 0.60 fiber volume ratio. The piezoceramic PZT5A was selected as the 
material for the piezoelectric layer(s) with properties provided by the manufacturer (ref. 35). This material is 
recommended for sensor applications by the manufacturer. The material properties are shown in Table 1. 

Eight ply symmetrically laminated beams with a single delamination and 2 piezoceramic layers, each 
bonded to the top and bottom surface of the composite beam, were considered in the numerical examples. 

Unless otherwise stated, a cross-ply layup [ p/0/ 90/0/90] s was mostly considered, where “p” denotes a 
piezoceramic layer. Laminations of the type (p/0/90/45/-45] s and [p/0/45/-45/0] s were also investigated. The 
length of the beam was equal to 5 in, its width was 1 in, and the thickness of each layer was 0.01 in. Results 
are presented for a cantilever beam clamped at the left end. 

The beams were assumed to be excited by a uniformly distributed load of 1 lb/in amplitude applied on the 
upper surface of the beam at a frequency of 80 Hz which is lower than the fundamental frequency. This low 
frequency was selected to avoid excitation of the higher modes that could complicate the interpretation of the 
results. The terminals were assumed closely spaced having infinitesimal length, such that they yield continuous 
voltage readings along the free surfaces. 


Effect on Dynamic Deflections 

The first phase of the analysis evaluates the effect of delamination on the transverse dynamic deflections of 
the sensory beam. In the case of a cantilever beam these effects are quite significant. As follows from figure 2, 
increases in the delamination length resulted in significant increases in the amplitude of vibration. This is 
associated with the reduction of stiffness in the delaminated beam. The predicted influence of local vibration 
modes in the delaminated regions on the global mode shape appears negligible both in the present and in the 
following examples. 
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The location of a delamination through the depth of the cantilever beam had a minimal effect on the motion 
(see fig. 3). The two cases compared in figure 3 include: a delamination located at the mid-plane of the beam 
(mid-plane case); and a delamination located near the upper surface (off-plane case) such that [p/0/90/0/90/90/0- 
//90/0/p], where // indicates the presence of the crack. Considering various axial locations of delaminations in 
figure 4, the axial location of a delamination has a relatively larger effect. As expected, the displacement 
amplitudes increase as the delamination approaches the clamped end of the beam, indicating that such 
delaminations result in more pronounced reductions of the stiffness. 

Sensory Response 

The following paragraphs present the predicted voltage output across the piezoelectric sensors. The outputs 
from the cantilever beam with centrally located delaminations of various lengths are shown in figures 5 to 7. In 
these figures, as well as in the following figures, the two curves represent the amplitudes of the resultant AC 
voltage across the lower and upper piezoelectric layer (referred to as “lower” and “upper” layers) respectively. 
Apparently the presence of the delamination can be detected either from the resultant discontinuities in the 
voltage output along each surface, or from the resultant difference in the voltage amplitudes between the sensors 
in the upper and lower surfaces (see figs. 5 to 7). More importantly the size and axial location of the delamina- 
tion are sharply outlined. In practice, such an accurate detection of a delamination will be limited by the finite 
number and size of sensors/electrodes, however, the spacing of the sensors is a design problem which should be 
solved in conjunction the allowable minimum length of a detectable delamination. 

The predicted effect of an axially shifted delamination on the voltage outputs of a cantilever beam is shown 
in figures 8 to 9. In both cases, voltage readings clearly indicate the presence, axial location and size of the 
delamination. The resultant voltage change is lower when the delamination is near the ffee-end, as opposed to a 
delamination near the root. This happens because the axial strain in the piezoelectric layers also decreases 
towards the free end. 

The voltage output of the piezoelectric sensors is also sensitive to the through-the-thickness location of 
delamination (see figs. 10 and 11). Figure 10 corresponds to a delamination located near the top surface of the 
beam, while figure 1 1 to a delamination near the bottom surface. The differences between the figures (compare 
also to figure 5 where the delamination lies on the mid-plane) indicates that the proposed technique may also 
reveal the lateral location of a delamination, in addition to detecting the presence and axial location. 

The laminate layup did not affect significantly the conclusions obtained in the previous examples. This was 
verified for a central delamination in cantilever beams with laminations [p/0/90/45/-45] s and [p/0/45/-45/0] s 
respectively. Figure 12 compares the predicted voltage output of a [p/0/90/0/90] s and a [p/0/45/-45/0] s cantilever 
beam with 20 percent central delaminations at their mid-plane. Although there are some differences in the 
voltage output, the detection of delamination is not affected. 

The proposed technique is not also affected by the supporting conditions of the beam. This is clearly demon- 
strated in figure 13, where the sensory response of a [p/0/90/0/90] s simply supported beam is shown. The beam 
has a 20 percent central delamination at the mid-plane. As was expected the different supporting conditions 
result in different voltage distributions, yet, the presence, size and axial location of delamination may be 
detected by the discontinuities in the voltage distributions along each surface. Additional case studies were 
conducted, similar to the ones presented for the sensory cantilever beam case, and indicated that the simply- 
supporting conditions din not affect the effectiveness and capability of the proposed technique to detect the 
presence, size, and axial/lateral location of a delamination. 
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SUMMARY AM) CONCLUSIONS 


A technique for detecting delaminations in vibrating composite beams was proposed, and the associated 
mechanics were developed. The proposed scheme is based on monitoring the changes in the voltage output 
signature of distributed piezoelectric sensors bonded to or embedded into the vibrating structure. The vibrations 
may be induced either externally (mechanically or acoustically) or by using piezoelectric actuators. 

The method was applied to detect delaminations in composite beams. Based on the obtained results, the 
following conclusions are summarized. The predicted response of the piezoelectric sensors was very sensitive to 
the presence of a delamination. It was observed that in the case of continuous piezoelectric layers with 
infinitesimally spaced terminals bonded to the upper and lower surfaces, abrupt axial discontinuities in the AC 
voltage output across the sensors indicate clearly the crack tips. Moreover, resultant voltage differences between 
the sensors in upper and lower surfaces may be used to directly detect the presence of delamination. Hence, the 
proposed concept can detect the presence of a delamination and predict its size mid location both axially and 
through-the-thickness of a composite structure. The effectiveness of the detection was not affected by the type 
of boundary conditions. This was illustrated by analyzing simply supported and cantilever beams. 

In closing, the reported research has demonstrated the feasibility of the proposed technique. The method 
described herein can be successfully extended to detect delaminations in composite plates and shells. Future 
work will address the effectiveness of the method when limited numbers of sensors are implemented, and the 
experimental verification of the proposed technique. 
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APPENDIX 


Elements of Matrix [ay] and Vector {fj} for a Composite Beam with Single Delamination. 


The right end of the beam is free to move in the axial direction while the left end is completely restrained. 
Therefore, it can be shown that/^ = V 4 = 0 and N 2 = - N 3 . The system of equations (19) is reduced to 21 

equations. The vectorjcj includes the following elements: 



p f 1 z'-' 1 f 1 s-,2 pi q2 

L l C 2 C 3 C 4 C 5 Cj C 2 C 3 C 4 


cf c\ cl cl cl cl cl cl cl cl n 2 j 


I. Simply supported beam subjected to a mechanical load 


Matrix [ay]: The elements a ij are written in the form a(i,j). Elements that are not listed are equal to zero. 

1 . i- 1 

fl(l,2)= a(l,3)= a(l,4)= 1 

2. i- 2 

a(2,2)= - a(2,3) = - a{ 2,4)= D l 

3. i- 3 

a(3,l)=a(3,3)=-a(3,4)=B 1 a(3,5)= 1 

4. i= 4 

a(4,16)= t 4 a(4,17)= c 4 a(4,18)= e 4 a(4,19)= e w 

6. i= 6 

a(6,l)= - B n c x 

5. i= 5 a(6,5)= - 1 

a(5,16)= D 4 f 4 a(5,17)= D 4 c 4 a(5,18)= - D 4 e 4 a(5,19)= - D 4 e m 

a(6,9)= - B 2 e 2 Q 


a(7,l)= - Bj 3 Cj 
a(7,5)= - 1 

a (7,14)= ^3^30 


o(7,2)= ^13^1 
a(7,ll)=-B 3 c 3 

a(7,15)= 1 


<*(7,3) = ^i3 e i 
<*(7,12) = - S 3 f 3 

<*(7,21) = - li 
A 3 


a(7,4)= B j 3 e 10 
a(7,13)=5 3 e 3 


8. i= 8 

a(8,l)=r, <*(8,2) = c, <*(8,3)=^ 

a(8,6) = - t 2 <*(8,7) = - c 2 <*(8,8) = - e 2 


<*(8,4) = e 10 
<3(8,9) = - e 20 
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9. i-- 9 

a( 9,1)= /j o( 9,2)= q o(9,3) = e x a( 9,4)= e 10 

0(9,11)= - > 3 o(9,12)= - c 3 o(9,13)= - e 3 o(9,14)= - e 30 

10. i= 10 

o(10,l)= JjCj o(10,2)= - Sjfj o(10,3) = s l e l o(10,4)= - jjg 10 

o(10,6)= - s 2 c 2 a(10,7)= s 2 t 2 o(10,8)= - s 2 e 2 a(10,9)= s 2 e 20 

11. i- 11 

o(ll,l)=5 1 c 1 0(11,2)= - o(ll,3)= Sjgj a(ll,4)= - jje 10 

o(ll,ll)= - s 3 c 3 o(ll,12)= s 3 t 3 o(ll,13)= - s 3 e 3 a(ll,14)= s 3 e 30 

12. i -- 12 

o(12,6)= B 2 c 5 a(12,7)= - B 2 t 5 o(12,8)= B 2 e 5 o(12,9)= - B 2 e 50 

o(12,10)= 1 a(12,16)= - B n c 6 o(12,17)= B n t e o(12,18)= - B n e 6 

I, 

o(12,19)= B X2 e 6 o o( 12,20)= - 1 o(12,21)= _ 

13. i- 13 

a(13,ll)= B 3 c 7 o(13,12)= - 5 3 r 7 o(13,13)= B 3 e 7 a(13,14)= - B 3 e 10 

o(13,15)= 1 a(13,16)= - B l3 c 6 a(13,17)= B x3 t 6 o(13,18)= - B X3 e 6 

I , 

o(13,19)= B^e^ a( 13,20)= - 1 a(13,21)= - _ 

A 3 

14. i= 14 

o(14,6)= t 5 o(14,7)= o 5 o(14,8)= e 5 o(14,9= e 50 

a(14,16)= - t 6 a(14,17)= - c 6 o(14,18)= - e 6 a( 14,19)= - 

15. i- 15 

a(15,ll)= r 7 a(15,12)= c 7 a(15,13)= e 7 a(15,14)= e 1Q 

o(15,16)= - t 6 o(15,17)= - c 6 a(15,18)= - e 6 o(15,19)= - 

16. i- 16 

a(16,6)= - s 2 c 5 a(16,7)= s 2 t 5 a(16,8)= - s 2 e 5 a(16,9)= s 2 e 50 

a(16,16)= SjC 6 o(16,17)= - s x t 6 o(16,18)= SjC 6 a(16,19)= - 

17. i = 17 

a(17,ll)= - s 3 c 2 a(17,12)= s 3 / 7 a(17,13)= - s 3 e 7 a(17,14)= ^ 3 e 70 

o(17,16)= 5 jC 6 a(17,17)= - s x t§ a(17,18)= s x e 6 a(17,19)= - s x e^ 

18. i= 18 

a(18,l)= - DjSj 3 Cj a(18,2)= D x s x t x a(18,3)= D x s x e x o(18,4)= - D x s x e xo 

o(18,6)= D 2 s 2 c 2 o(18,7)= - D 2 s 2 t 2 a(18,8)= - D 2 s 2 e 2 a(18,9)= b 2 s 2 e 2G 

o(18,ll)= D 3 s 3 c 3 a(18,12)= - D 3 s 3 t 3 a(18,13)= - D 3 s 3 e 3 a(18,14)= D 3 s 3 e 3Q 

19. i= 19 

a(19,6)= D 2 s 2 c 5 o(19,7)= - D 2 s 2 t 5 o(19,8)= - D 2 s 2 e 3 o( 19,9)= b 2 s 2 e 50 

a(19,ll)= D 3 s 3 c 7 a(19,12)= - b 3 s 3 t 7 o(19,13)= - b 3 s 3 e 7 o(19,14)= D 3 s 3 e 1Q 

a(19,16)= - b x sfc 6 o(19,17)= D x s x t 6 o(19,18)= D x sfe 6 a(\9,l9)- - 
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20. i- 20 


a(20,l) = D l Syt l <2(20,2) -- D x s x c x 

a( 20,6)= - D 2 s^t 2 <2(20,7)= - D 2 s 2 c 2 

a(20,ll) = - D 3 5 3 2 2 3 <z(20, 12) = - D 3 s 3 c 3 

B-y B’J 

a(20,21) = — - — - z 2 + ^3 
A 2 A 3 


<a(20,3)= - 
a(20,8)= D 2 s 2 e 2 
<2(20,13= D 3 s 3 e 3 


21 . /= 21 

a(21,6)= - D 2 s 2 2 r 5 

a(21,ll) = - D 3 5 3 2 r 7 

a(21,16)= 

a(21, 21)= a(20,21) 
In the above equations, 

B n s n 


B n- 


In addition. 


a(21,7)= - D 2 s 2 c 5 
a ( 21 , 12 )= - D 3 s 3 c-i 
<2(21,17)= b { s x c 6 


B ' 


Dn" »n- 




S. 




a(2 1,8)= D 2 s 2 e 5 
a(21,13)= D 3 s 3 e-j 

a(21,18)= - Dj5 2 e 6 


<2(20, 4 ) = - iy 2 e 10 
<2(20,9)= b 2 s 2 e 20 
<2(20,14)= D 3 s 3 e 30 


a( 21,9)= D 2 5 2 e 50 
<2(21,14)= D 3 s 3 e 2XX 
<2(21,19)= - b x s x e(fi 


c r - cos S k l^ 
S *1 


t r - sin s k \ q 


e r --e 


e rO z e 




A correspondence between subscripts r, k and q is established in the following Table 


r 

i 

2 

3 

4 

5 

6 

7 

k 

i 

2 

3 

1 

2 

1 

3 

q 

i 

1 

1 

! q = L 

2 

2 

2 


B. Vector {fj } 

Elements that are not listed are equal to zero. 

ft = f 4 = - Fj 
f 8 = ‘ f l4 = F 2 ' F i 
f 9 = * f 15 = F 3 _1 F 1 

II. Cantilever beam clamped at the left end 

The second and the fourth lines of the matrix [ay] are changed as follows: 
2 = 2 

<2(2,1) = <2(2,3) = - a(2, 4) = 1 
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ft ■ fs ■ • ^- b iK 

A \ 



i - 4 

a( 4,16) = - i>jSj 3 c 4 a{4,Yl) = DjSj 3 ^ 

a(4,18) = D l s^e 4 a( 4,19) = - D l s^e 4Q 

The element f 4 of the vector {fj} is equal to zero, contrary to the case of a simply supported beam. 
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TABLE L — MATERIAL PROPERTIES 


[In material coordinates] 


Property 

T300/934 

PZT-5A 
(Morgan Matroc, 
(ref. 35)) 

E E ,,(GPa) 

Vy 0.275 
E~ 2 2 (GPa) 

G 66 (GPa) 
d,, (mAO 
e S 33 (farad/m) 
p (kg/m 3 )14957700 
Ply Thickness(mm) 

150 

0.275 

7.3 

3.2 

0 

0 

1495 

0.2544 

69 

0.31 

69 

26.3 

-171 10~ 9 

15.045 10" 9 
7700 
0.254 



Delaminations 



A B D C A 

* «r — t* -r — — " 


(b) 

Figure 1 . — Delaminated sensory beam, (a) Typical configuration; (b) Typical beam segments. 
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Figure 2. — Effect of delamination length on the dynamic displace- 
ment of a cantilever beam. Centrally located delamination. 



Normalized Length (x/L) 

Figure 4. — Effect of the axial location of a delamination on the 
dynamic displacement. Delamination length = 20%. Coordi- 
nate of delamination center is x^L Left: x^L = 0.3; Right: 
x^L = 0.7; Center: x^L = 0.5. 



Figure 6.— Detection of a centrally located delamination in a 
cantilever beam. Delamination length = 40%. 



Figure 3. — Effect of through-the-thickness location of a delamin- 
ation on the dynamic displacement. Delamination length = 20%; 
lamination: [p/(0/90>2/90/0//90/0/p]. 



Figure 5. — Detection of a centrally located delamination in a 
cantilever beam. Delamination length = 20%. 



Normalized Length (x/L) 

Figure 7.— Detection of a centrally located delamination in a 
cantilever beam. Delamination length = 60%. 
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Normalized Length (x/L) 

Figure 8. — Detection of an axially shifted delamination in a canti- 
lever beam. Delamination length = 20%. Delaminatlon center 
is located at 30% span from the clamped end. 



Normalized Length (x/L) 


Figure 1 0. — Detection of a delamination located near the top 
surface of a cantilever beam. Delamination length is 20%; 
off-plane delamination: [p/(0/90)2/90/0//90/0/p]. 



Normalized Length (x/L) 

Figure 12. — Effect of laminate configuration on delamination 
detection. Laminate configurations: [p/0/45/-45/0], and 
[p/0/90/0/90] s . 20% delamination length centrally located at 
mid-plane. 



Figure 9. — Detection of an asymmetric delamination in a canti- 
lever beam. Delamination length = 20%. Delamination center 
is located at 70% span from the clamped end. 



Normalized Length (x/L) 

Figure 11 . — Detection of a delamination located near the bottom 
surface of a cantilever beam. Delamination length is 20%; 
off-plane delamination: [p/0/90//0/90/(90/0);j/p]. 



Figure 1 3. — Effect of simply supporting conditions on delamination 
detection. Surface of a cantilever beam. [p/Q/90/0/90] s beam 
with 20% delamination length centrally located at mid-plane. 
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